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Abstract 



In this paper, we consider near cloaking for the full Maxwell equations. We extend the 
method of where the quasi-static limit case and the Helmholtz equation are considered, 

to electromagnetic scattering problems. We construct very effective near cloaking structures 
for the electromagnetic scattering problem at a fixed frequency. These new structures are, 
[ before using the transformation optics, layered structures and are designed so that their first 

scattering coefficients vanish. Inside the cloaking region, any target has near-zero scattering 
cross section for a band of frequencies. We analytically show that our new construction 
, significantly enhances the cloaking effect for the full Maxwell equations. 
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^ ; 1 Introduction 

QQ ' The cloaking problem is to make a target invisible from far-field electromagnetic measurements 

! [ISl HOI 1221 [SI HH m] • Many schemes for cloaking are under active current investigation. These 

' include exterior cloaking in which the cloaking region is outside the cloaking device [25l [1] [24j [TTl 

P\J I Unilllj active cloaking [15], and interior cloaking, which is the focus of our study. 

In interior cloaking, the difficulty is to construct electromagnetic material parameter distribu- 
tions of a cloaking structure such that any target placed inside the structure is undetectable to 
waves. One approach is to use transformation optics [351 [131 [301 [331 [IS]- It takes advantage of 
the fact that the equations governing electromagnetism have transformation laws under change of 
KJi • variables. This allows one to design structures that steel waves around a hidden region, returning 

j_j ' them to their original path on the far side. The change of variables based cloaking method uses 

a singular transformation to boost the material properties so that it makes a cloaking region look 
like a point to outside measurements. However, this transformation induces the singularity of 
material constants in the transversal direction (also in the tangential direction in two dimensions), 
which causes difhculty both in the theory and applications. To overcome this weakness, so called 
'near cloaking' is naturally considered, which is a regularization or an approximation of singular 
cloaking. In [19] . instead of the singular transformation, the authors use a regular one to push 
forward the material constant in the conductivity equation describing the quasi-static limit of elec- 
tromagnetism, in which a small ball is blown up to the cloaking region. In [18| . this regularization 
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point of view is adopted for the Helmholtz equation. See also ^2^1 HH] • More recently, Bao and Liu 
[5] considered near cloaking for the full Maxwell equations. They derived sharp estimates for the 
boundary effect due to a small inclusion with an arbitrary material parameters enclosed by a thin 
high-conducting layer. Their results show that the near cloaking scheme can be applied to cloak 
targets from electromagnetic boundary measurements. 

In [SI [3] it is shown that the near cloaking, from measurements of the Dirichlct-to-Neumann 
map for the conductivity equation and of the scattering cross section for the Helmholtz equation, 
can be drastically enhanced by using multi-layered structures. The structures are designed so 
that their generalized polarization tensors (GPTs) or scattering coefficients vanish (up to a certain 
order). GPTs are building blocks of the far-field behavior of solutions in the quasi-static limits 
(conductivity equations) and the scattering coefficients are 'Fourier coefficients' of the scattering 
amplitude. The multi-layered structures combined with the usual change of variables (transfor- 
mation optics) greatly reduce the visibility of an object. This fact is also confirmed by numerical 
experiments [3]. 

The purpose of this paper is to extend the results of [Sj [6] to Maxwell's equation and show 
that the near cloaking from cross section scattering measurements at a fixed frequency can be 
enhanced by using layered structures together with the change of variables. Again the layered 
structures are designed so that their first scattering coefficients vanish. It is also shown that inside 
the cloaking region, any target has near-zero scattering cross section for a band of (low) frequencies. 
We analytically show that our new construction significantly enhances the near cloaking effect for 
the full Maxwell equations. It is worth mentioning that even if the basic scheme of this work is 
parallel to that of [6 , the analysis is much more complicated due to the vectorial nature of the 
Maxwell equations. 

The paper is organized as follows. In Section 2, we recall some fundamental results on the 
scattering problem for the full Maxwell equations. In Section 3, we introduce the scattering 
coefficients of an electromagnetic inclusion and prove that the scattering coefiicients are basically 
the spherical harmonic expansion coefficients of the far-field pattern. Section 4 is devoted to the 
construction of layered structures with vanishing scattering coefficients. We also present some 
numerical examples of the scattering coefficient vanishing structures. In Section 5, we show that 
the near cloaking is enhanced if a scattering coefficient vanishing structure is used. 

2 Multipole solutions to the Maxwell equations 

In this section, we recall a few fundamental results related to electromagnetic scattering, which 
will be essential in the sequel. 

Consider the time-dependent Maxwell equations 



where u) is the frequency. The couple (E, H) is a solution to the harmonic Maxwell equations 




where fx is the magnetic permeability and e is the electric permittivity. 

In the time-harmonic regime, we look for the electromagnetic fields of the form 




V X E = icj^H, 

V X H = -iweE. 



(2.1) 



We say that (E, H) is radiating if it satisfies the Silver-Miiller radiation condition: 

lim |x|(77IH X i - ViE) = 0, 
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where x — x/|x|. In the sequel, we set k = uj^jlji, which is cahed the wave number. 

For m = — n, . . . ,n and n = 1,2,..., set Y^^ to be the spherical harmonics defined on the unit 
sphere S. For a wave number A: > 0, the following function 

^;„,™(fc;x) = /iW(fc|x|)r™(x) (2.2) 

satisfies the Helmholtz equation At; + fc^u = in M"^ \ {0} and the Sommerfeld radiation condition: 

lim |x|(^p^(fc;x) - ifcw„,„(fc;x)) = 0. 

|x|->-oo CIXl 

Here, /in ■* is the spherical Hankel function of the first kind and order n which satisfies the Som- 
merfeld radiation condition. Similarly, z;„,m(x) is defined as 

<m(A;;x)=j„(A;|x|)r™(x), (2.3) 

where j„ is the spherical Bessel function of the first kind. The function Vn,m satisfies the Helmholtz 
equation in all R"^. 

In the same manner, we can make solutions to the Maxwell system with the vector version of 
spherical harmonics. Define the vector spherical harmonics as 

^ VsFr(*) and V„,„=xxU„,„, (2.4) 



for m — —n^ . . . ,n and n = 1, 2, . . . . Here, x S S* and V5 denotes the surface gradient on the unit 
sphere S. The vector spherical harmonics defined in (|2.4I) form a complete orthogonal basis for 
L|,(S'), where L'^{S) = {u €E (i^(S'))^ | f • u = 0} and v is the outward unit normal to S. 

Multiplying the vector spherical harmonics to the Hankel function, we can make the so-called 
multipole solutions to the Maxwell system. To make the analysis simple, we separate the solutions 
into transverse electric, (E • x) = 0, and transverse magnetic, (H • x) =0. Define the exterior 
transverse electric multipoles to p.ip as 



E^,^.(fc;x) - ~v/^(n + l)/iW(fc|x|)V„,,„(x), 



H^,i(fc;x) = Vx - v/n(n+l)/iW(fc|x|)V„,„,(x) 

LUfl \ 



(2.5) 



and the exterior transverse magnetic multipoles as 



E™(fc;x) = ^Vx( - v/n(n + l)/iW(fc|x|)V„,™(x)), 
H™(fc;x) = -v/n(n + l);iW(fc|x|)V„,„(x). 

The exterior electric and magnetic multipole satisfies the radiation condition. By the same way, 
we define the interior multipoles (EJ^,HJ_^) and (Ej_^^,H^^^) with hn^ replaced by jn, i.e., 



E^.^(fc;x) = -v/n(n + l)j«(fc|x|)V„,„(i), 

H^.t(fc;x) = --VxE^,^(fc;x), ^'-^^ 

UJfl 

and _ 

H™(fc;x) = -Vn(n + l)j«(fc|x|)V„,^(x), 

E^*^(fc;x) = ^VxH™(fc;x). ^^'^^ 

ute 

We will sometimes omit the wave number k in the notation of the multipoles. 
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Note that we have 

V X E^,^(fc;x) = ^/^j|+IlH„(fc|x|)U„,„(x) + !^i!J±ll/,(i)(fco|x|)K„™(x)x, (2.9) 

V X E^,t(fc;x) - ^"f'|^^V „(fc|x|)u„,„,(x) + !^i!i±llj(i)(fco|x|)rr(x)x, (2.10) 

|x| |x| 

where = h';^\t) + t (/li^^)' (t) and J'„(t) j„(t) + tj'^{t). 

The solutions to the Maxweh system can be represented as separated variable sums of the 
multipole solutions, see [27l Section 5.3]. With multipole solutions and the Hclmholtz solutions in 
(|2.2p and (|2.3p . it is also possible to expand the fundamental solution to the Helmholtz operator. 
For fc > 0, the fundamental solution to the Helmholtz operator (A + k'^) in M.^ is 



rfe(x) 



47r X 



(2.11) 



Let p be a fixed vector in R^. For |x| > |y|, the following addition formula holds (see [26l Section 
9.3.3]): 



Ffe(x-y)p = -^^^- J2 E™(fc;x)EjM(fc;y).p 
— ' n(n + 1) a ^ — ' 

n— 1 ^ ^ m— — n 

ik 



+ T.;j^) E E^,i(fc;x)Ej,^(fc;y).p 

n—1 m— — n 

OO 71 



k 



(2.12) 



Ti— 1 m— — n 



with Vn^m and "Dn.m being defined by (|2.2p and ()2.3p . 

Plane wave solutions to the Maxwell equations have the expansion using the multipole solutions 
as well (see [TT]). The incoming wave E*(x) = z/c(q x p) x qe^'^'^ '', where q g is the direction of 
propagation and the vector p S is the direction of polarization, is expressed as 



E^(x)=zfc^. 



where c = (q x p) x q. 



(Vp,,(q) . c)Ejl(x) - — (Up,,(q) • c)Ej,f (x) 



(2.13) 



3 Scattering coefficients of an inclusion 

Let Z? be a bounded domain in R"^ with C^'" boundary for some a > 0, and let (eo, /io) be the pair 
of electromagnetic parameters (permittivity and permeability) of M'^ \ and (ei, /ii) be that of D. 
We assume that eo,ei,/io, and fii are positive constants. Then the permittivity and permeability 
distributions are given by 

e = eox(R'\^)+eix(I?) and ^ = /iox(R' \ ^D) + /iix(i?), 

where x denotes the characteristic function. In the sequel, we set k = cj^ei/ii and fco = w^egMo- 
For a given solution (E*, H*) to the Maxwell equations 

r V X E^ = iw^off in R^ 
\ V X H* -iweoE^ in R^, 
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let (E, H) be the solution to the following Maxwell equations: 

V X E = iujfiU in R^, 

V X H = -icueE in R^, (3.1) 
(E — E*, H — H') satisfies the Silver-Miiller radiation condition. 

We emphasize that along the interface dD, the following transmission condition holds: 

[i/ X E] = [i/ X H] = 0. (3.2) 
Here, [v x E] denotes the jump of x E along dD, namely, 

[z/ X E] = (i/ X E)|+^ - (j/ X E)|g^. 
Let Vac' denote the surface divergence. We introduce the fimction space 
TH{div,dD) :=|m e L^dD) : Vqd ue L^idD)^ 
equipped with the norm 

\\u\\TH(div,dD) = \\u\\l2(0D) + ||Va_D ' u\\L2f^Qjjy 

For a density cp £ TH{div, dD), we define the single layer potential associated with the fundamental 
solutions Ffe given in (|2.11l) by 



cS^[<p](x) := / r,(x-y)<p(y)da(y), xG 
JdD 



For a scalar density contained in L^{dD), the single layer potential is defined by the same way. 
We also define boundary integral operators: 



£^[<^](x) -.^{ux (es^^] + V5|j[VaD • v])) (x), 
X^[<p](x) :=p.v. / i.(x)x(Vxx(Ffe(x-y)<^(y)))dfT(y), ^eOD. 

JdD ^ ' 

In the same way, we define S^^ , C^^ , and associated with F/j^ instead of Ffc. Then the solution 
to (|3.ip can be represented as the following: 



E(x) 



and 



E*(x) + X [¥j](x) + V X V X [i/;](x), x e \ D, 
AtiV X 5^[v3](x) + V X V X 5^[V'](x), x e D, 



H(x) = (V X E) (x), X e N 



(3.3) 



where the pair (^p, ip) e TiJ (div, dD) x Ti? (div, dD) is the unique solution to 



rk r>ko 



rk r-ko 
'-D '-D 

2fj,i 2/io/ /ii /io 







■ E'' X 1/ ■ 








iojH* X u 





dD 



(3.4) 

The invertibility of the system of equations p. 41) on TH{div, dD) x riJ(div, 9_D) was proved in 
|32) . Moreover, there exists a constant C = C{e,iJ.,ijj) such that 



l¥'llTff(div,az5) + ll'0llTH(div.aD) < C(||E* x vWrHidi^MD) + l|H' 



(3.5) 



From f|2.12p (with /cq in the place of k) and (|3.3p it follows that, for sufficiently large |x|, 

oo . , n 

(E-E^)(x) = ^--^1^ ^ (a„,„E^,^(fco;x) + /3„,„E^,^(A:o; x)) , (3.6) 



^ n(n+ 1) ^ ^ 

n— 1 m——n 

-zc^eoMo / E^(fco; y) • ^(y) + kl f WKiko; y) • ^^(y), 

JdD JdD 



where 



I3„,m = -iujeof^o / Ej|^(fco; y) • ^(y) - 4 / ^n^i^o; y) • ip{y). 

JdD JdD 

Definition 1. Let (Vp.f , ) be the solution to Jg.^p iw/ien E* = Ep^(fco; y) and H* = Hj_^(fco; y), 
anrf {'f™,i>™) when E' = E™{ko;y) and W = U™{ko]y). The scattering coefficients 
(wT^'^^f VFJ^^'T.*^ -..wT^^'T,^ \,WT^''Tf^ a associated with the permittivity and the perme- 

\ [n,m}{p,q}^ (n,m)(p,g)' {n,m)(p,q)^ {n,m)(p.q) J r ti r 

ability distributions e,/i and the frequency lj (or k,kQ, D) is defined to be 

^mL) = / ^^(''o; y) • Vp,f (y) My) + kl f WKiko; y) • ^Pl^iy) da{y), 

J BD J dD 

TE,TM _ . ,, f ^TMIi. ...\ .,TAn,.\ , 1,2 f ^tk li, . ..\ .i.TM 



Wln%tL) - -'^'ol^o / K'^iko; y) • (y) da{y) + ki / Ejfjfco; y) • V^,f (y) da{y), 

J dD J dD 



W^K^Kp,,) = -*^^oA^o / Ej,^,(fco;y)-<^(y)d^(y)-^^e^ / ^™{k,;y) ■ rp::;{y) da{y), 

JdD JdD 

= -^^^oA^o / E^(fco; y) • v™{y) da{y) - Lo^el [ E^(fco; y) • V^f (y) da(y). 

dD J dD 

As we see it now, the scattering coefficients appear naturally in the expansion of the scattering 
amplitude. We first obtain the following estimates of the scattering coefficients. 

Lemma 3.1. There exists a constant C depending on {e,fi,uj) such that 



rTrTE,TE r 1 
^(n.m)(p,g)[e'A^'H 



(jn+p 



for all n,m,p,qe N. The same estimates hold for W^J^^J,,,, ^™n^ip,q)' '^'^'^ 

Proof. Let (v','0) be the solution to dSJ]) with E*(y) = Ej^(fco;y) and H* = -^V x W. Recah 
that the spherical Bessel function jp behaves as 



l-3---(2p+l) V \p, 
uniformly on compact subsets of M. Using Stirling's formula p\ — \/27rp(p/e)^(l + o(l)), we have 

f CPtP\ 

JpW^Of— asp^oo, (3.8) 

uniformly on compact subset of R with a constant C independent of p. Thus we have 

flip 

1 1 ■pi 1 1 _|_ ||"fj*|| <' 

II \\TH(div,dD) II \\TH{div,dD) — pp 



for some constant C . It then follows from (13. 5p that 

II II II /II — 

W'^WL^idD) + ll''^'llL2(a_D) - 

for another constant C. So we get p.7p from the definition of the scattering coefficients. □ 
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Suppose that the incoming wave is of the form 

oo p 



(3.9) 



p=l q=-p 

for some constants a^.g and bp ^. Then the solution (<p, -0) to (|3.4p is given by 

oo p 

TM 



p=l q=-p 



^ = E E y^P^l'^p'q + ^P.'J'/'J 



TA/ 



By (I3.6|) and Definition [1] the sohition E to p.ip can be represented as 



oo . , n 

(E-E^)(x) = ^--^^^ ^ (a„,,„E^,^(fco;x) + /3„,„E™(A:o; x)) , |x| ^ oo, (3.10) 

n— 1 ^ m——n 



where 



oo p 



^n,m — E/ E/ ('^P 
p=l g=-p 
oo p 

E E 



y^TE,TE ^ y^TE,TM 

^ (n.m){p^q) {n,m)(p.q) 



TTrTM,TE , „.TM,TM 

^^{n,m){p,q) (n,m)(p,q) 



(3.11) 



p=l q=-p 

Using (|3.10p . (13. lip and the behavior of the spherical Bessel functions, we can estimate the 
far-field pattern of the scattered wave (E — E'). The far-field pattern (also called the scattering 
amplitude) Aoo[e,Ai,w] is defined by 



gifco|x| 

E(x) - E'(x) = Aoo[e,A^,^](x) -|- o(|x|"^) as |x| oo. 



/co|x 



(3.12) 



Since the spherical Bessel function hn"^ behaves like 



as t — ^ oo, 



as t — > oo. 



one can easily see by using ()2.9p that 



E^^t(A:o;x) ^ ^e-^- ( - ^^^(^rTT)) V„,™(i) 



El^(fco;x) 



fco|x 



fco|x| V eo 
The following result holds. 



-e ' 2 



as X — OO, 



as X oo. 



Proposition 3.2. // E' is given by I13.9\) . then the corresponding scattering amplitude can be 
expanded as 



oo '^'^/c j 

Aoo[e,^,cj](x) = V' ° (a„.„V„.„(x) -h/3„,™./— U„,„i(x) 

^1 V^('^ + l) ™t^„^ V eo 

where an,m and Pn.m are defined by h3.11]) . 



(3.13) 
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We emphasize that since {Vn,m, ^n,m} forms an orthogonal basis of L'^{S), the conversion of 
the far-field to the near field is achieved via formula (|3.10p . 

We now consider the case where the incident wave is given by a plane wave e^'^ ^'c with 
|k| = fco and k • c = 0. It follows from (I2TT3)) that 



Vp,,(k) • c)Ej,f (fco;x) - . (Up,,(k) • c)Ej*^(fco;x 



where k = k/fco G -S*, and therefore 
4niP 



(Vp,g(k)-c) and bp^, 



(Up,,(k)-c) 



Hence, the scattering amplitude, which we denote by Aoo[e,M, w](c, fc; a;), is given by p.l3p with 
the coefhcients a„.m and Pn.m 



EE 



AniP 



p=l q=-p 

Pn,m — 2^ I—, — V 

p=l ,= -p ^JV{V + 1) 



(Vp,,(k).c)W-, 



TE,TE 
(n,m)(p,q) 



^ lu li \ \'sxrTE,TM 

(Up,g(k)-c)VK(„,„)(p,^) 



1 



(Up.,(k).c)W-J*X^) 



(3.14) 

These formulas tell us that the scattering coefficients appear in the expansion of the scattering 
amplitude. 

We now investigate the low frequency behavior of the scattering coefficients. Let r(x) := 
— l/(47r|x|) denote the fundamental solution corresponding to the case k = 0, and Aio the associ- 
ated boundary integral operator: 

Xz,[¥>](x) :=p.v. / i.(x)x(Vxx(r(x-y)<^(y)))d(7(y), ^eTH{dw,dD). 

JdD ^ ' 

Analogously to p.4p . one can prove that there is a unique solution (y5(°' , e TH{div, dD) x 

TH{div, dD) to the following equations: 



Ml + Mo 
2(Aii - Mo) 




I + Md 



(ei - eo) 





(1 + (0 



■I + Mi 







' E * X ■ 













dD 



.2(ei - eo) 

(3.15) 

In fact, since dD is C^'", Aio is compact and we may apply the Fredholm alternative to prove 
unique solvability of above equation. Moreover, we have 



'\\TH(di^,dD) +i^W '\\TH(div,dD) < C(||E' X uWrHidiv.dD} + ||H' X v\\TH{div,dD)), (3.16) 

with a constant C ~ C(e, m)- 

Let M be a small positive number and consider the boundary integral equation (|3.4p with fc, ko, 
and uj replaced by pk, pko, and puj, respectively. Then, we have (see jl4l ) 



M'j^-MD = Oip^), M'j^° - Md = Oip') 



and 
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Since 



(^ + 3. 
V2/Z1 2/zo 



-M 



D 



Mo 



ei + eo 



I + {ei-eo)MD + 0{p'') 



if we express the solution to p.4p as (ipTip) ■— [v'' ^ P^'4>^) ^ then it satisfies 



A + 0(p) 



iH' X u 



dD 



where A is the 2-by-2 matrix appeared in the left-hand side of p.lSp . From the invertibility of A, 
it follows that there are constants po and C — C(e, p,, cj) independent of p as long a.s p < po such 
that 



\\v''\\TH{div,dD} + P^\\'4>''\\TH{div,dD} < C(||E' X l^\\TH{div,dD) 

Lemma 3.3. There exists po such that, for all p < po, 

(jn+p 



TH(div,dD)) 



rrTTE,TE r 1 

^(n,m)(p,g)[e'A^'H 



< 



n+p+l 



(3.17) 



(3.18) 



for all n,m,p,q G where the constant C depends on {e,p,uj) but is independent of p. The same 
estimate holds for W^^ ^^^^^^^y W^(n,m)(p,9)' a'^'^ W^(„,m)(p,g) ■ 

Proof. Let (yj.t/j) be the solution to dSH) with E*(y) = 'e'^^^ {pko;y) and H' = x E*. 



Then, from 



it follows that 

IIe'^^II 



\TH(div,dD) 

where C is independent of p, and hence 



, ■ 11 CP 

\ \\TH(diM,dD) - pP ^ ' 



< 



CP 



for p < Pa for some pQ. So we get p.lSp from the definition of the scattering coefficients in 
Definition m □ 



4 S-vanishing structures 

The purpose of this section is to construct multilayered structures whose scattering coefficients 
vanish, which we call S-vanishing structures. The multi-layered structure is defined as follows: For 
positive numbers ri, . . . , r^+i with 2 — ri > r2 > ■ ■ ■ f'L+i = 1, let 

A, {x : rj+i < |x| < r,}, j = l,...,L, Ao:=R^\A;, Al+i{= D) := {x : \x\ < I}, 

and 

r, = {|x|-r,}, j = l,...,i + l. 

Let (/ij, (.j) be the pair of permeability and permittivity parameters of Aj for j = Set 
^0 = 1 and £0 = 1- We then define 

L+l L+l 

^=Z!^J^(^J") ^^'^ e = I]ejX(^j), (4.1) 

j=0 j=0 
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which are permeability and perniittivity distributions of the layered structure. 

The scattering coefficients ( WT^^^f s^WT'^'^.f^ ^,Wj''^'Z'^ ',,WT^Z'^ are defined as be- 

\ {n,m){p,q}^ {n,'m){p^q)^ (n,m)(p,g)' {n,m}{p,q} J 

fore, namely, if E* given as in p.9p . the scattered field E — E* can be expanded as (|3.10p and 
p. lip . The transmission condition on each interface Tj is given by 



[i X E] = [i X H] = 0. 
We assume that the core A^+i is perfectly conducting (PEC), namely, 

E X IV = on r^+i = dA^+i. 



(4.2) 



(4.3) 



Thanks to the symmetry of the layered (radial) structure, the scattering coefficients are much 
simpler than the general case. In fact, if the incident field is given by E' — E^^, then the solution 
E to (|311) takes the form 



E(x) 



ElmW + ajE/,_t„(x), xeAj, j=0,...,L, 



(4.4) 



with flQ = 1. From (12.91) and (I2.10p . the interface condition (|4.2p amounts to 



1 



Jn{kj-ir.j) 



1 



T-Ln{kj-irj) 



a J -I 



(4.5) 



where Unit) = h'it\t) + 1 {h!:;}^^' {t) and Jn{t) = jn{t) + tj'Jt), and the PEC boundary condition 

(4.6) 



on r 



V„(fcL) hi'\kL) 




aL 




0" 







aL 








L + l IS 



Since the matrices appeared in (j4.5l) are invertible, one can see that there are aj and aj, j = 
0, 1, . . .L satisfying (|4.5I) and (14. 6p . Similarly, one can see that if the incident field is given by 
E' = E^*;^(x), then the solution E takes the form 



E(x) = h,Y.l%{^) + 5,E^;^(x), X e A„ J = 0, 1, L 



T^TM 



(4.7) 



for some constants bj and bj {bo = 1). One can sec now from 
coefficients satisfy 



and (|4.7p that the scattering 



w, 



/ \r \ — W, \, for all (171,71) and (v.q), 

{n^m){p,q) {n,m){p^q) ^ ' ^ KfiHj-i 

(„,m)(p,9) = ^(n,m)(p,g) = if (m, fl) ^ {p, q) , 



and, since (|4.4p and (|4.7p hold independently of to, we have 



W, 



W, 



TE,TE 
(n,0)(n,0) 
TM.TM 



w, 



(n,0)(n,0) 



W, 



TE,TE 
(n,m) (n,m) ' 
TM.TM 



(n,m) (n,m) 



for — n < 717 < n. 



Moreover, if we write 



W'^^ — W'^^ and l/F^*^ — 
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then we have 



wTB tn{n+l) j.^ in{n+l) 
= ao and = 

fco fco 



60. 



(4.8) 



Suppose now that E^{^ is the incident field and the solution E is given by 
E(x) = 5,E^|(x) + ajE^_^(x), x e A„ j = 0, . . . , i, 

with ao = 1, where the coefficients Oj's and Oj's are determined by (|4.5p and 
(1431) that 



. We have from 



1 



1 



for j = 1, 



where 



xn 

We then have from 



L. Substituting these relations into (|4.6p yields 



(4.9) 



Pn.l Pn,2 


- 1 



Mi 





■'Hn{kjrj 



-;ili)(fc,r,)' 



jn(fcL) 









'Jn{kjrj) jn{kjrj) 

fj,j 



T-Lnikj-lTj) 



(4.10) 



m(n + 1) 



Mt(n + 1) p^f 



-ao 



(4.11) 



Similarly, for W^f*^, we look for another solution E of the form 

E(x)=6,E^*^(x) + 6,E^*^(x), xe^,, 3 = 0,...,L, 
with bo ~ 1. The transmission conditions become 



Jn{kjrj) T~Ln{kjrj) 



jn{kj 



1 



h^n\k,r,) 



1 



£,_1 



Hn{kj-irj) 



and the PEC boundary condition on the inner most layer is 



.^ + 1, 



JnikL) n^ikL) 







'bL 




0" 




h. 








From (|4.12p and (|4.13|) . we obtain 



bo 
bo 



(4.12) 



(4.13) 



(4.14) 
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where 



■ TM TM 
Pn,l Pn,2 









1 



^— J«(fcj-irj) -^Hnikj^iTj) 
in{kj-irj) hl^\kj.irj) 



(4.15) 



From the definition of and (|4.14l) , we have 



,TM 



w. 



TE 



in{n + 1) bo m(n + 1) p„.i 



fco 



'^O Pr, 



(4.16) 



It sliould be emphasized tliat p„f ^ and p„ 2 7^ 0- f^*^^' if 2 ~ tlien (|4.9p can be 
fulfilled with aq = and oq = 1. This means that there exists (/i, e) on K.'^ \ D such that the 
following problem has a solution: 

V X E = iw/iH in \ 

V X H = -iojeE in \ 



(x X E) 







on dD, 



E(x) = E^,f(x) for|x|>2. 
Applying the following Green's theorem on 57 = {x | 1 < |x| < R}, 

/ (E • AF + curlE • curlF + divE divF)dx 
Jn 

= / (i^ X E • curlF + • E divF)dcr(x) 
Jdn 



with F = E„ q(x) and the PEC boundary condition on {x| — 1}, we have 



x|=_R, 



(i^ X E) • Hrfcr(x) = iko / (|H|^ - lE^rfx 



In particular, the left-hand side is real-valued. Hence, 



|H X jy- Epdcr(x) 



|x| = fi 



(|H X + |Ep - 25R((i/ X E) • H)dcr(x) 
[ (|Hxj.|2 + |E|2)da(x). 

J\x\=R 



From the radiation condition, the left-hand side goes to zero as i? — oo, and it contradicts the 
behavior of the hankel functions. One can show that p^^^ ^ in a similar way. □ 

To construct the S-vanishing structure at a fixed frequency to we look for (/.t, e) such that 



TMr 



for some N. More ambitiously we may look for a structure (/i, e) for a fixed w such that 
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for all 1 < n < and p < po for some po- Such a structure may not exist. So instead we look for 
a structure such that 

Wr[f^,e,pu] = o(p2^+i), Wr'[p,e,P^] = (4.17) 

for all 1 < n < A'^ and p < po for some po. We call such a structure a S-vanishing structure of 
order N at low frequencies. In the following section, we expand the scattering coefficients for low 
frequencies and derive conditions for the magnetic permeability and the electric permittivity to be 
a S-vanishing structure. 

Suppose that {p, e) is an S-vanishing structure of order N at low frequencies. Let the incident 
wave E* be given by a plane wave e^^^'^c with |k| = fcg and k-c = 0. From (|3.14p . the corresponding 
scattering amplitude, Aoo[p, e, pw](c, fc; x), is given by (|3.13p with the following an,m and /3„,m- 

47r7" 

-.{V^,^{k)-c)W^^[p,e,pu;], 



Jn{n -f 1) 

47r7" 1 

/3„,™ = , , (U„,,„(k) • c)W^™[A.,e,H- 

+ 1) «w^o 

Applying (|3.18p and (|4.17p . we have 

Aoo[Ai,e,pH(c,k;x)=o(p2^+i) (4.18) 

uniformly in (k,x) if p < Po- Thus using such a structure the visibility of scattering amplitude is 
greatly reduced. 

4.1 Asymptotic expansion of the scattering coefficients 

The spherical Bessel functions of the first and second kinds have the series expansions: 



XI 2'ni- 3- • -(271 + 2^ + 1) 
and 



(2n)! ^ (_i)'t2/-n-i 



2"n! ^ 2n\{-2n + l)(-2n + 3) • • • (-2n + 2? - 1) ' 

So, using the notation of double factorials, which is defined by 

f n • (n - 2) . . . 3 • 1 if n > is odd, 
n!! := < n • (n — 2) . . . 4 • 2 if n > is even, 
[ 1 ifn = -l,0, 

we have 

■?-^(0- (^^^^),, (l + o(0) for^«l, (4.19) 

and 

?/„(0 = -((2n- l)!!)t-"+i(l + o(0) fort<l. (4.20) 
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We now compute P^^[e,/i,t] for small t. For n>l 



(2n+ 1)!! 











. /i,(2n+l)!! * ^''^^ ^ 



(2n+ 1)!! 



r + o(r) 



-«Q(n) 



L Mi-i(2n+l)!! ' > M.,--ife-ir-,)"+^' + / 



where Zj = ^JeJWj and Q{n) = (2n — 1)!!. We then have 



(2n + l)!! ^ ^ ^ ^2+' 



n 



(2n+l)..„^. 

n n + l 2ri + ] 

' ((2n + l)!!) 
Similarly, for the transverse magnetic case, we have 

T->TM r .1 



"^^""^("+^) (1 - J^)t-+^ (1 + (n + 1 + ^) (1 + 



n 

3 = 1 



(2n + l)!! + ^ ^ ^2+' 



Q(n)«5'_i e, , ,,,, , (<3(n))^n 



((2n + l)!!)2 I' e,_J* + 



2n + l 



e,_i 



r^"-! (1 + 0(1)) 



,n + l 



Q{n)z 
{2n + l)\\z';+l 



n + 1 + (1 + 0(1)) 

{'"E „„J „T_E 



Using the behavior of spherical bessel functions for small arguments, we see that p^f and p^f 
admit the following expansions: 

(N-n \ 
5:/„7(M,e)i^'+o(t")j 



and 



/N-n 



\ 1=0 J 

Similarly, "p^^ and "p^^ have the following expansion: 



and 



fN-n 



I iV —11 \ 

for t = pijo and some functions fn,f^9nf^f™^ 9™ independent of t. 



(4.21) 

(4.22) 

(4.23) 
(4.24) 
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Lemma 4.1. For any pair of e), we have 
and 



(4.25) 



(4.26) 



Proof. Assume that there exists a pair of (/^, e) such that g^^lji^e) — 0. Since [/i,e, pw] = 
o(p~"~^), the solution given by (|4.4p with ao = 1 and ao = satisfies 

V X X - p^w^eE = in \ 

V • E = in R3 \ 15, 
(i/ X E)|_^ = o(p-("+i)) on 91), 

I E(x) = /ili)(pfco|x|)y„^o(x) for |x| > 2. 

Let V(x) = linip^o p"^^E(x). Using ()4.20p we know that the Umit V satisfies 



V X ( -V X V 

V • V = 







in R''' \ D, 



in R3 \ D, 

(iyxy)|^=0 on^D, 
[ y(x) = -{{2n - l)!!)K,o(A) for |x| > 2. 

Since V^,o(x) = 0(|a;|~"'^), we get V(x) = by Green's formula, which is a contradiction. Thus 
5n (a^j ^) 7^ 0- ^ similar way, (|4.26p can be proved. □ 



From Lemma |4. 11 we have the following theorem. 
Proposition 4.2. We have 



N-n 



N-n 



rjrTMr ^ ,i ,2n+l TjrTMr t,21 . f,2N+l\ 



where t ^ puj and the coefficients W^f[fj,,e] and W^^f'^ [fx, e] are independent oft. 
Hence, if we have {fi, e) such that 

Wlf [/i, e] = W™ [p, £] = 0, for aU 1 < n < TV, < ; < (TV _ n). 



(4.27) 



(/X, e) satisfies (|4.17p . in other words, it is a S-vanishing structure of order N at low frequencies. 
It is quite challenging to construct (/i, e) analytically satisfying (|4.27p . In the next section we get 
some numerical examples of such structures. 



4.2 Numerical examples 

In this section we provide numerical examples of S-vanishing structures of order TV at low frequen- 
cies based on (j4.27p . To do this, the gradient descent method for the suitable energy functional is 
used, as used in j5j and [S] to compute the enhanced near-cloaking structures for the conductivity 
problem and the Helmholtz problem. As in [51, we symbolically compute the scattering coefficients. 
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In the place of spherical Bessel functions and spherical Hankel functions, we put its low frequency 
asymptotic expansions in (|4.10l) and (I4.15|) . and symbolically compute W^^ and W^^^ to have 
Wlf[^^,e]^nd W™[^^.e]. ^ 

The following example is a S-vanishing structure of order N — 2 made of 6 multilayers. The 
radii of the concentric disks are Vj = 2 — for j = 1, . . . , 7. From Proposition 14. 2[ the nonzero 
leading terms of W^^[fi, e, t] and e,t] up to are 

• [t^,t^] terms in W^f^[^,e,i], i.e., VF^^ff , W'l^f , 

• [^3,^5] terms in W^^' e,t], i.e., W™ ,W™ , 

• [t^] term in W^^[y.,e,t\, i.e., W^^ , 

• [t^] term in e,t], i.e 
Consider the mapping 



"2,0 ■ 



1.0 , , W^i n" , Wff, Wli, Wm 



(4.28) 



where, /i = (/xi, . . . , /xg) and e = (ei, . . . , e%). We look for (/x, e) which has the right-hand side of 
(|4.28p as small as possible. Since (|4.28p is a nonlinear equation, we solve it iteratively. Initially, 
we wet /i. = /i.'-"-' and e — e*^"^ We iteratively modify (/J*-*-', e*-'-*) 



where A\ is the pseudoinverse of 



(i+l) „(»+l)lT _ 



(4.29) 



A,; 



1 J • 



9(/X,£) 



(M,e) = (M<'',eW) 



and 



(M,^) = (M<'':^<'') 



Example 1. Figure UTT] and Figure W?^ show computational results of 6- layers S-vanishing structure 
of order N ^2. We set r = (2,^,...,|), /x^^) = (3,6,3,6,3,6) and e(°) = (3,6,3,6,3,6) and 
modify them following (|4.29p with the constraints that and e belongs to the interval between 0.1 
and 10. The obtained material parameters are ^l = (0.1000,1.1113,0.2977,2.0436,0.1000,1.8260) 
and e = (0.4356,1.1461,0.2899,1.8199,0.1000,3.1233), respectively. Differently from the no-layer 
structure with PEC condition at |x| = 1, the obtained multilayer structure has the nearly zero 
coefficients oiW^^[^,e,t\ and W^'^''[n,e,t\ up to . 



5 Enhancement of near cloaking 

We make a cloaking structure based on the following lemma. 

Lemma 5.1. Let F he a dijjeomorphism of M?' onto R? such that i^(x) is identity for |x| large 
enough. If (E, H) is a solution to 

' V X E = 2w/iH in R?, 

< V X H = -iueE m E^, (5.1) 

(E — E*, H — H*) is radiating, 
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Figure 4.1: This figure shows the graph of the material parameters and the corresponding coef- 
ficients in W^^[iJ,,e,t] and W™[^,€,t] up to t^. The first row is of the no-layer case, and the 
second row is of 6- layers S-vanishing structure of order N = 2 which is explained in Example 1. In 
the third column, the y-axis shows {W^^ ,Wlf ,Wl^ ,Wlf ,Wl^^ ,W^/) from the left to the 
right. 




12345 12345 

n n 



Figure 4.2: This figure shows the graph of W^^[ii,€,t] and W^^[ii,e,t] for various values of t. 
The first row is of the no-layer case, and the second row if of 6-layers S-vanishing structure of order 
N = 2 which is explained in Example 1. The values of W^^ and are much smaller in the 

S-vanishing structure than in the no-layer structure. 

then (E,H) defined by (E(y),H(y)) =(E(F-i(y)), H(F-i(y))) satisfies 

V X E = iuj{F^ij:)ii.^ m M^, 

V X H = -iw(F*e)E in M?, 
(EE*, H - H*) ^5 radiating, 
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where (E^(y),ff(y)) ^{^\F-\y)),W{F-\y))) , 

, DFhL)ahL)DF^{-K) , , , DF(yi)e(yi)DF^ (yi) 

dctpF(x)) ' (^--)^^)= detpF(x)) ' 

it)it/i X = i^^^(y) BKC? is the Jacobian matrix of F. 
Hence, 

A[^,e,uj] = A[F^,fi, F^e,uj]. 

To compute the scattering amplitude which corresponds to the material parameters before the 
transformation, we consider the following scaling function, for small parameter p, 

4'i(x) = ix, xeM^. 

p 

Then we have the following relation between the scattering amplitudes which correspond to two 
sets of differently scaled material parameters and frequency: 

^oo fio ^ i^eo ^ i^uj ^ Aoo[p,e, puj]. (5.2) 
p p J 

To see this, consider (E, H) which satisfies 

(V X E) (x) = iw(/xo *i)(x)H(x) forxeM3\B^, 

(V X H)(x) = -iw(eo4i)(x)E(x) for x e \ B^, 

X X E(x) = on dBp, 
(E - E*, H - is radiating, 

with the incident wave E'(x) = e^'^'^'c and = x E' with k • c = and Ikl = kn. Here Br, 

is the ball of radius p centered at the origin. Set y ~ -x and define 



E(y),H(y)) :=((E o *^l)(y), (H o vl/^i)(y)) =((E o v]/^) (y), (H o (y] 
(y),ff (y)) :=((E^ ° *p)(y), (ff o vl/^)(y)^ 



E* 



Vy X E ) (y) = ipc^M(y)H(y) for y e \ Bi 



and 
Then, we have 



^Vy X Hj (y) - -zpt^e(y)E(y) for yeW\Bi, 

y^xE(y)^0^ ondBi, 
[ (E-E\H-H*) is radiating 

Remind that the scattered wave can be represented using the scattering amplitude as follows: 

gifco|x| 



(E-E^)(x) 



fco|x| 



(c, k; x) as |x| -> oo. 



and 



(E - E') (y) - - — p-r^°° ^' '^l as |y I ^ cx). 

KoP|y| 



Since the left-hand sides of the previous equations are coincide, we have (|5.2I) . 

Suppose that (/x, e) is a S- vanishing structure of order N at low frequencies as in Section 21 
From (|i?TS)) and Ij^^ . we have 



(c,k;x)=o(p2^+i) 



(5.3) 
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Then, we define the difFeomorphism Fp as 

; for Ixl > 2, 



f — 

V2(l 



3 - 4p 1 , A X 

'x| — for2p<|x|<2, 



^(l-p) 4(1 -p)' 7|x| 

+ forp<|x|<2p, 

— for |x| < p. 
P 



We then get from (|5.3p and Lemma [5. II the main result of this paper. 

Theorem 5.2. // (/i, e) is a S-vanishing structure of order N at low frequencies, then there exists 
po such that 

Aoo [(Fp)4^o*i),(Fp),(eo*i), J (c,k;i) = o(p2^+i), 
for all p < po, uniformly in (k, x). 

Remark that the cloaking structure ({Fp)^,{p,o ^t^i), (i^p),(e o vj/ 1 )) in Theorem 15.21 satisfies the 
PEC boundary condition on |x| = 1. 

6 Conclusion 

We have shown near-cloaking examples for the Maxwell equation. We have designed a cloaking 
device that achieves enhanced cloaking effect based on the method of O [3] to electromagnetic 
scattering problems. Any target placed inside the cloaking device has an approximately zero scat- 
tering amplitude. Such cloaking device is obtained by the blow up using the transformation optics 
of a multi-coated inclusion with PEC boundary condition. The cloaking device has anisotropic 
permittivity and permeability parameters. 
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